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K-matrix and Dalitz plot analysis from FOCUS 
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Dalitz analysis is a powerful tool for physics studies within and beyond the Standard Model. In the last decade 
it has helped to investigate the Heavy Flavor hadronic decay dynamics and is now being applied to extract 
angles of the CKM Unitarity triangle. To perform such sophisticate analyses we need to model the strong 
interaction effects. The FOCUS experiment has performed pilot studies in the charm sector through the K- 
matrix formalism. What has been learnt from charm will be beneficial for future accurate beauty measurements. 
Experience and results from FOCUS are presented and discussed. 
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1. Introduction 

Over the last years we have seen a resurrection of 
Dalitz plot analyses in modern Heavy Flavor experi- 
ments. This analysis tool, first applied in the charm 
sector, has more recently become a standard tech- 
nique, used for sophisticated studies and searches for 
new physics in the beauty sector. Paradigmatic exam- 
ples are B — > pit and B — > D^K^ for the extraction 
of the a and 7 angles of the Unitarity Triangle. In- 
deed, the road to go from the detected final states 
to the intermediate resonances can be rather insid- 
ious and complications arise in both decays. More 
precisely, the extraction of a in B — > pit means, op- 
eratively, selecting and filtering the desired interme- 
diate states among all the possible (tttt)tt combina- 
tions, e.g. an, /o(980)7r etc. The extraction of 7 
in B — > D^K^*) requires, in turn, modeling the D 
amplitudes. This poses the problem of how to deal 
with strong-dynamics effects, in particular those re- 
garding the scalar mesons. The mr and Kn S-wave 
are characterized by broad, overlapping states: uni- 
tarity is not explicitly guaranteed by a simple sum of 
Breit-Wigner functions. In addition, independently of 
the nature of the a, it is not a simple Breit-Wigner. 
The /o(980) is a Flatte-like function, and its lineshape 
parametrization needs a precise determination of KK 
and 7T7T couplings. Recent analyses of CP violation 
in the B — > DK channel from the beauty factories 
have used the Cabibbo-favored mode K s tt + tt~ , which 
is common to both D° and D°. A set of 16 two- 
body resonances had to be introduced to describe the 
(Ktt)tt and K s (ttw) states in the D° amplitude: two 
ad hoc resonances were required to reproduce the ex- 
cess of events in the tttt spectrum, one at the low-mass 
threshold, the other at 1.1 GeV 2 . Masses and widths 
of the two states, named <j\ and 02, were fitted to the 
data themselves and found to be M ai — 484 ± 9 MeV, 
r ffl = 383±14and M a . 2 = 1014±7MeV, = 88±13 
in BaBar [1] and M ai = 519 ±6 MeV, T ai = 454 ± 12 
M CT2 = 1050±8MeV, r CT2 = 101 ±7 in Belle [2]. These 
scalars were invoked with no reference to those found 
in other processes, in particular scattering data, and 
with no assumption as to the correctness of the physics 
the model embodies. This procedure of "effectively" 



fitting data invites a word of caution on estimating 
the systematics of these measurements. A question 
then naturally arises: in the era of precise measure- 
ments, do we know sufficiently well how to deal with 
strong-dynamics effects in the analyses? 

We have faced parametrization problems in the 
FOCUS experiment and learnt that many difficul- 
ties are already known and studied in different fields 
of physics, such as nuclear and intermediate-energy 
physics, where broad, multi-channel, overlapping res- 
onances are treated in the K-matrix formalism. The 
effort we have had to make mainly consisted in build- 
ing a bridge of knowledge and language to reach the 
high-energy community; our pioneering work in the 
charm sector might inspire future accurate studies in 
the beauty sector. 



2. The K-matrix and P-vector formalism 

A formalism for studying overlapping and many- 
channel resonances was proposed long ago and is 
based on the K-matrix [3, 4] parametrization. The 
K-matrix formalism provides a direct way of imposing 
the two-body unitarity constraint, which is not explic- 
itly guaranteed in the simple sum of Breit-Wigners, 
here referred to as the isobar model. Minor unitar- 
ity violations are expected for narrow, isolated reso- 
nances but more severe ones exist for broad, overlap- 
ping states. This is the real advantage of the K-matrix 
approach: it heavily simplifies the formalization of any 
scattering problem since the unitarity of the S matrix 
is automatically encoded. 

Originating in the context of two-body scattering, 
the formalism can be generalized to cover the case of 
production of resonances in more complex reactions 
[5], with the assumption that the two-body system 
in the final state is an isolated one and that the two 
particles do not simultaneously interact with the rest 
of the final state in the production process [4]. The 
validity of the assumed quasi two-body nature of the 
process of the K-matrix approach can only be verified 
by a direct comparison of the model predictions with 
data. In particular, the failure to reproduce the Dalitz 
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plot distribution could be an indication of the presence 
of relevant, neglected three-body effects. 



3. The FOCUS results 

3.1. The three pion analysis 

The FOCUS collaboration has implemented the K- 
matrix approach in the D s and D + — ► tt + tt^tt + anal- 
yses. It was the first application of this formalism in 
the charm sector. Results and details can be found in 
[6]. Here I only reprodece plots of the final results. In 
Fig. 1 and Fig. 2 the Dalitz-plot projections are shown 
for D s and D + into three pions. 



parameters was needed. The confidence levels of the 
final fits are 3.0 % and 7.7 % for the D s and D + re- 
spectively. The results were extremely encouraging 
since the same K-matrix description gave a coherent 
picture of both two-body scattering measurements in 
light-quark experiments as well as charm-meson de- 
cay. This result was not obvious beforehand. Fur- 
thermore, the same model was able to reproduce fea- 
tures of the D + — » 7r + 7r~7r + Dalitz plot that would 
otherwise require an ad hoc a resonance. The bet- 
ter treatment of the S'-wave contribution provided by 
the K-matrix model was able reproduce the low-mass 
7r + 7r~ structure of the D + Dalitz plot. This suggests 
that any cr-like object in the D decay should be con- 
sistent with the same cr-like object measured in -k + tt~ 
scattering. 




Low m 2 / s - mass projection (GeV 2 /c 4 ) High m 2 s + B - mass projection (GeV 2 /c 4 ) 



Figure 1: FOCUS Df Dahtz-plot projections with fit re- 
sults superimposed. The background shape under the signal 
is also shown. 




Low m 2 K + n - mass projection (GeV 2 /c 4 ) High 



mass projection (GeV 2 /c 4 ) 



Figure 2: FOCUS D Dahtz-plot projections with fit re- 
sults superimposed. The background shape under the signal 
is also shown. 

In Fig. 3 the FOCUS adaptive binning schemes for 
Df and D + are plotted. In this model [5], the pro- 
duction process, i.e, the D decay, can be viewed as 
consisting of an initial preparation of states, described 
by the P-vector, which then propagates according to 
(I — iKp)^ 1 into the final one. The K-matrix here 
is the scattering matrix and is used as fixed input in 
our analysis. Its form was inferred by the global fit 
to a rich set of data performed in [7]. It is interest- 
ing to note that this formalism, beside restoring the 
proper dynamical features of the resonances, allows 
for the inclusion in D decays of the knowledge coming 
from scattering experiments, i.e, an enormous amount 
of results and science. No re-tuning of the K-matrix 





Figure 3: D s and D + adaptive binning Dalitz-plots for the 
three pion FOCUS K-matrix fit. 

Further considerations and conclusions from the 
FOCUS three-pion analysis were limited by the sam- 
ple statistics, i.e. 1475 ± 50 and 1527 ± 51 events for 
D s and D + respectively. We considered mandatory 
to test the formalism at higher statistics. This was 
accomplished by the D + — > K~7t + tt + analysis. 



3.2. The D+ K 



The recent FOCUS study of the D+ -► K-tt+tt+ 
channel uses 53653 Dalitz-plot events with a signal 
fraction of ~ 97%, and represents the highest statis- 
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tics, most complete Dalitz plot analysis for this chan- 
nel. Invariant mass and Dalitz plot are shown in Fig. 4. 
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Figure 4: Fhe D + — » if ~7r + 7r + Dalitz plot (left) and mass 
distribution (right): signal and sideband regions are indi- 
cated in red and blue respectively. The sidebands are at 
±(6-8) a from the peak. 



the same K-matrix formalism. However, the contri- 
bution from the vector wave comes mainly from the 
if * (892) state, which is well separated from the higher 
mass if* (1410) and K* (1680), and the contribution 
from the spin-2 wave comes from if.* (1430) alone. 
Their contributions are limited to small percentages, 
and, as a hrst approximation, they can be reasonably 
described by a simple sum of Breit-Wigners. More 
precise results would require a better treatment of 
the overlapping if* (1410) and if* (1680) resonances 
as well. In accord with SU(3) expectations, the cou- 
pling of the if 7T system to if 77 is supposed to be sup- 
pressed. Indeed we find little evidence that it is re- 
quired. Thus Fx/2 i s actually a vector consisting of 
two components: the first accounting for the descrip- 
tion of the if 7r channel, the second of the if r}' chan- 
nel: in fitting D + — > if ~ir + ir + we need, of course, the 
(F 1 / 2 ) 1 element. Its form is 



Details of the analysis are in [8] . 

An additional complication in the Kit system comes 
from the presence in the S*-wave of the two isospin 
states, I = 1/2 and I = 3/2. Although only the 
I = 1/2 is dominated by resonances, both isospin 
components are involved in the decay of the D + me- 
son into if ~7T + 7r + . A model for the decay amplitudes 
of the two isospin states can be constructed from the 
2x2 K-matrix describing the 1 = 1/2 S-w&ve scat- 
tering in (-KV)i and (ifr/)2 (with the subscripts 1 
and 2, respectively, labelling these two channels), and 
the single-channel K-matrix describing the 1 = 3/2 
K~ir + — > if ~7r + scattering. 

The K-matrix form we use as input describes the S- 
wave if ~7r + — > if ~7r + scattering from the LASS ex- 
periment [9] for energy above 825 MeV and K~ir~ —* 
K~ir~ scattering from Estabrooks et al. [10]. The 
K-matrix form follows the extrapolation down to the 
if7r threshold for both I = 1/2 and 1 = 3/2 S- wave 
components by the dispersive analysis by Biittiker et 
al. [11], consistent with Chiral Perturbation Theory 
[12]. The complete form is given below in Eqs. (4-5) 
with the parameters listed in Table I [13]. 

The total -D-decay amplitude can be written as 

M = (i^W+^OO + E^ ^ B (^), (!) 

3 

where s = M 2 (Kir), (F 1 / 2 ) 1 and F 3 / 2 represent the 
I = 1/2 and 1 = 3/2 decay amplitudes in the Kir 
channel, j runs over vector and spin-2 tensor reso- 
nances 1 , and B{abc\r) are Breit-Wigner forms. The 
J > resonances should, in principle, be treated in 



(F 1/2 ) 1 = (I- l K 1/2 p)^ 1 (P 1/2 ) J , (2) 

where I is the identity matrix, if 1/2 is the K-matrix 
for the I = 1/2 S- wave scattering in if7r and if 77', p 
is the corresponding phase-space matrix for the two 
channels [4] and (Pi/ 2)] is the production vector in 
the channel j. 

The form for F3/2 is 

F 3 / 2 = (I-iK 3/2P )- 1 P 3/2 , (3) 

where if 3/2 is the single-channel scalar function de- 
scribing the I = 3/2 K~n + — > if ~7r + scattering, and 
P3/2 is the production function into if7r. 

Fitting of the real and imaginary parts of the 
if ~7r + — » if ~7r + LASS amplitude, shown in Fig. 5, 
and using the predictions of Chiral Perturbation The- 
ory to continue this to threshold, gives the K-matrix 
parameters in Table I. 



LASS LASS 




M (Kit) (GeV) M (Kit) (GeV) 

Figure 5: Real and imaginary if~7r + — > K~ir + ampli- 
tudes from the LASS experiment and their K-matrix fit 
results. 



1 Higher spin resonances have been tried in the fit with both 
formalisms but found to be statistically insignificant. 



The 1=1/2 K-matrix is a single-pole, two-channel 
matrix whose elements are given in Eq. (4) . 
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Kn = ('-J-i) + C 110 + C m 5 + C 112 S 2 ) 

\ ^norm / \ s l s / 

K 22 = ( S —^-) ( 9 -^J2l + C 220 + C 221 s + C 222 ~s 2 ) 

\ s norm / \ «1 - S J 

K 12 = ( ^i) ( ^ + c 120 + c 121 s + c 122 sA , 

\ ^norm / \ ^1 5 _/ 

(4) 

where the factor of s norm = m 2 K + m 2 is conve- 
niently introduced to make the individual terms in 
the above expression dimcnsionlcss. g\ and 172 are 
the real couplings of the s\ pole to the first and the 
second channel respectively. s Q i = 0.23 GeV 2 is the 
position of the Adler zero in the 1=1/2 ChPT elas- 
tic scattering amplitude 2 . Cm, C 22 i and C\ 2 i for 
i = 0,1,2 are the three coefficients of a second or- 
der polynomial for the diagonal and off-diagonal el- 
ements of the symmetric K-matrix. Polynomials are 
expanded around § = s/s norm — 1. This form gener- 
ates an S-matrix pole, which is conventionally quoted 
in the complex energy plane as -E = M — iT/2 = 
1.408 - iO.110 GeV. Any more distant pole than 
Kg (1430) is not reliably determined as this simple 
K-matrix expression does not have the required an- 
alyticity properties. Nevertheless, it is an accurate 
description for real values of the energy, where scat- 
tering takes place. Numerical values of the terms in 
Eq. (4) are reported in Table I. 

The 7 = 3/2 K-matrix is given in Eq. (5). Its form 
is derived from a simultaneous fit to LASS data [9] 
and to 7r — > K~n~ scattering data [10]. It is a 
non-resonant, single- channel scalar function. 

K V2 = ( S -^-l \ (D uo + D ni s + D 112 S 2 ) . (5) 

\ Snorm / 

In Eq. (5) s 3 = 0.27 GeV 2 is the Adler zero position 
in the 7 = 3/2 ChPT elastic scattering and the values 
of the polynomial coefficients are -Dno = —0.22147, 
D m = 0.026637, and D 112 = -0.00092057 [13]. 

When moving from scattering processes to D- 
decays, the production P-vector has to be introduced. 
While the K-matrix is real, P-vectors are in general 
complex reflecting the fact that the initial coupling 
D + — ► (K~ ^ + )^t P ectator ncc d not be real. The P- 
vector has to have the same poles as the K-matrix, 
so that these cancel in the physical decay amplitude. 
Their functional forms are: 



2 Chiral symmetry breaking demands an Adler zero in the 
elastic 5-wave amplitudes in the unphysical region. ChPT at 
next-to-leading order fixes these positions Sqi [11, 12]. 



(A/ 2 )i = + (cio + c n s + c 12 s 2 )e^ (6) 



(^1/2)2 = + (c 20 + c 21 s + c 22 s 2 )e^ (7) 



^3/2 - (c 3 o + c 31 s + c 32 s 2 )e i73 . (8) 

(3e ie is the complex coupling to the pole in the 'ini- 
tial' production process, 171 and g 2 are the couplings as 
given by Table I. The Ktt mass squared s c = 2 GeV 2 
corresponds to the center of the Dalitz plot. It is con- 
venient to choose this as the value of s about which 
the polynomials of Eqs. (6-8) are expanded, by defin- 
ing s = s — s c . The polynomial terms in each channel 
are chosen to have a common phase 7j to limit the 
number of free parameters in the fit and avoid uncon- 
trolled interference among the physical background 
terms. Thus, the coefficients of the second order poly- 
nomial, aj, are real. Coefficients and phases of the P- 
vectors, except g\ and g 2 , are the only free parameters 
of the fit determining the scalar components. 

Free parameters for vectors and tensors are am- 
plitudes and phases (aj and Si). Kir scattering de- 
termines the parameters of the K-matrix elements 
and these are fixed inputs to this D decay analysis. 
Table II reports our K-matrix fit results. It shows 
quadratic terms in (-Pi/ 2 )i are significant in fitting 
data, while in both {P\/ 2 ) 2 and P 3 / 2 constants are 
sufficient. 

The J > states required by the fit are listed in 
Table III. 

The S'-wave component accounts for the dominant 
portion of the decay (83.23 ± 1.50)%. A significant 
fraction, 13.61 ± 0.98%, comes, as expected, from 
i^*(892); smaller contributions come from two vec- 
tors if* (1410) and if* (1680) and from the tensor 
if* (1430). It is conventional to quote fit fractions for 
each component and this is what we do. Care should 
be taken in interpreting some of these since strong 
interference can occur. This is particularly apparent 
between contributions in the same-spin partial wave. 
While the total S'-wave fraction is a sensitive measure 
of its contribution to the Dalitz plot, the separate fit 
fractions for I = 1/2 and 1 = 3/2 must be treated 
with care. The broad I = 1/2 S-w&ve component 
inevitably interferes strongly with the slowly varying 
1 = 3/2 ,5- wave, as seen for instance in [14]. Fit re- 
sults on the projections are re shown in Fig. 6. The 
corresponding adaptive binning scheme is at the top 
of Fig. 7 . 

The fit x 2 /d.o.f is 1.27 corresponding to a confi- 
dence level of 1.2%. If the 1 = 3/2 component is 
removed from the fit, the x 2 /d.o.f worsens to 1.54, 
corresponding to a confidence level of 10~ 5 . 
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Table I Values of parameters for the 1 = 1/2 K-matrix. 



pole (GeV 2 ) 


coupling (GeV) 


Cm 


Cl2i 


C22i 


si = 1.7919 


gi = 0.31072 
g 2 = -0.02323 


Ciw = 0.79299 
Cm = -0.15099 
C112 = 0.00811 


C120 = 0.15040 
C121 = -0.038266 
C122 = 0.0022596 


C220 = 0.17054 
C221 = -0.0219 
C222 = 0.00085655 



Table II S'-wave parameters from the K-matrix fit to the FOCUS D + — * K~tt + tt + data. The first error is statistic, 
the second error is systematic from the experiment, and the third is systematic induced by model input parameters for 
higher resonances. Coefficients are for the unnormalized S-wave. 



coefficient 


phase (deg) 


P = 3.389 ± 0.152 ± 0.002 ± 0.068 
cio = 1.655 ± 0.156 ± 0.010 ± 0.101 
en = 0.780 ± 0.096 ± 0.003 ± 0.090 
C12 = -0.954 ± 0.058 ± 0.0015 ± 0.025 
c 20 = 17.182 ± 1.036 ± 0.023 ± 0.362 
c 30 = 0.734 ± 0.080 ± 0.005 ± 0.030 


6 = 286 ± 4 ± 0.3 ± 3.0 

71 = 304 ± 6 ± 0.4 ± 5.8 

72 = 126 ±3±0.1±1.2 

73 = 211 ±10 ±0.7 ±7.8 


Total S-wave fit fraction = 83.23 ± 1.50 ± 0.04 ± 0.07 % 
Isospm 1/2 fraction = 207.25 ± 25.45 ± 1.81 ± 12.23 % 
Isospm 3/2 fraction = 40.50 ± 9.63 ± 0.55 ± 3.15 % 



Table III Fit fractions, phases, and coefficients for the J > components from the K-matrix fit to the FOCUS D + — > 
K~ir + ir + data. The first error is statistic, the second error is systematic from the experiment, and the third error is 
systematic induced by model input parameters for higher resonances. 



component 


fit fraction (%) 


phase Sj (deg) 


coefficient 


K*(892)tt+ 


13.61 ±0.98 
± 0.01 ±0.30 


(fixed) 


1 (fixed) 


iT(1680)7r+ 


1.90 ±0.63 


1±7 


0.373 ± 0.067 




± 0.009 ± 0.43 


±0.1 ±6 


± 0.009 ± 0.047 


if 2 *(1430)7r+ 


0.39 ±0.09 


296 ± 7 


0.169 ±0.017 




± 0.004 ± 0.05 


± 0.3 ±1 


± 0.010 ±0.012 


K*(UW)n+ 


0.48 ±0.21 


293 ± 17 


0.188 ±0.041 




± 0.012 ±0.17 


± 0.4 ±7 


± 0.002 ± 0.030 



These results can be compared with those obtained 
in the effective isobar model, which can serve as the 
standard for fit quality. Projections are shown in 
Fig. 8 and the adaptive binning scheme at the bot- 
tom of Fig. 7. 

Two ad hoc scalar resonances are required, of mass 
856 ± 17 and 1461 ± 4 and width 464 ± 28 and 177 ± 8 
McV/c 2 respectively. A detailed discussion of the re- 
sults and the systematics can be found in [8]. The 
results of the K-matrix fit showed that a consistent 
representation with scattering is possible, the global 
fit quality being indeed good. However, it deteriorates 



at higher Kir mass. This is not surprising since our 
K-matrix treatment only includes two channels Kir 
and Krj' . While we have reliable information on the 
former channel, we have relatively poor constraints 
on the latter. This means that as we consider Kir 
masses far above Krj threshold, these inadequacies 
in the description of the Krj' channel become increas- 
ingly important. This is expected to become worse as 
yet further inelastic channels open up. Consequently, 
improvements could be made by using a number of 
Z?-decay chains with Ktt final state interactions and 
inputting all these in one combined analysis in which 
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Adaptive binning and % contributions for dp to kpp (data) 



► mass projection (GeV 2 /c 4 ) 



m K - mass projection (GeV /(f) 



Figure 6: The Dalitz plot projections with the K-matrix 
fit superimposed. The background shape under the signal 
is also shown. 



;f"t - 202.114 
Free Par 16 
DOF = 159 

= 1.18E-02 
(Vdof = 1.271 
= 20 
i = 175 



1>X 2 <2 
2>X*<3 




several inelastic channels are included in the K-matrix 
formalism. In the present single D + — ► K~tt + tt + 
channel, adding further inelastic modes would be just 
adding free unconstrained parameters for which there 
is little justification. It is interesting to note that 
the adaptive binning scheme shows that both the K- 
matrix and the isobar fit are not able to reproduce 
data well in the region at 2 GeV 2 , in the vicinity of 
the K rj' threshold. It is also the energy domain where 
higher spin states live. Vector and tensor fit param- 
eters in the two models are in very good agreement: 
we do not exclude the possibility that a better treat- 
ment of these amplitudes could improve the \ 2 - Some 
isolated spots of high \ 2 could be caused by an im- 
perfect modeling of the efficiency as they are in the 
same regions in both fits. 

A feature of the K-matrix amplitude analysis is that 
it allows an indirect phase measurement of the sep- 
arate isospin components: it is this phase variation 
with isospin 1—1/2 which should be compared with 
the same 1—1/2 LASS phase, extrapolated from 825 
GeV down to threshold according to Chiral Perturba- 
tion Theory. This is done in the right plot of Fig. 9. In 
this model [5] the P-vector allows for a phase variation 
accounting for the interaction with the third particle 
in the process of resonance formation. It so happens 
that the Dalitz fit gives a nearly constant production 
phase. The two phases in Fig. 9b) have the same 
behaviour up to ~ 1.1 GeV. However, approaching 
Krj' threshold, effects of inelasticity and differing fi- 
nal state interactions start to appear. 

The difference between the phases in Fig. 9a) is due 
to the 1 = 3/2 component. 

These results are consistent with Kir scattering 
data, and consequently with Watson's theorem pre- 
dictions for two-body Ktt interactions in the low Ktt 
mass region, up to ~ 1.1 GeV, where elastic processes 
dominate. This means that possible three-body in- 
teraction effects, not accounted for in the K-matrix 
parametrization, play a marginal role. 

Our results for the total S- wave are in general agree- 
ment with those from the E791 analysis, in which the 



Adaptive binning and % contributions tor dp to kpp (data) 



X toe = 184 . 150 
Free PaF 16 
DOF = 157 

- 6.83E-02 
/dof - 1.173 
(fur = 20 
= 173 




Figure 7: The adaptive binning schemes corresponding to 
the K-matrix (top) and isobar (bottomt) fits. 



S-wave modulus and phase were determined in each 
Ktt slice [15], [16]. What does this analysis contribute 
to the discussion of the existence and parameters of 
the k! We know from analysis [17] of the LASS data 
(which in K~tt + scattering only start at 825 MeV) 
there is no pole, the k(900), in its energy range. How- 
ever, below 800 MeV, deep in the complex plane, there 
is very likely such a state. Its precise location re- 
quires a more sophisticated analytic continuation onto 
the unphysical sheet than the K-matrix representation 
provided here. This is because of the need to approach 
close to the crossed channel cut, which is not correctly 
represented for a robust analytic continuation. How- 
ever, our K-matrix representation fits along the real 
energy axis inputs on scattering data and Chiral Per- 
turbation Theory in close agreement with those used 
in the analysis by Descotes-Genon and Moussallam 
[18] that locates the K with a mass of (658 ± 13) MeV 
and a width of (557±24) MeV by careful continuation. 
These pole parameters are quite different from those 



Proceedings of the CHARM 2007 Workshop, Ithaca, NY, August 5-8, 2007 



7 




► mass projection (GeV 2 /c 4 ) 



m K - ,+ mass projection (GeV /c ) 



Figure 8: Dalitz plot projections with our isobar fit super- 
imposed. The background shape under the signal is also 
shown. 



closer to threshold than the older scattering results. 
This serves as a valuable check from experiment [20] 
of the inputs to the analyses of [18] and [19] based 
largely on theoretical considerations. Dalitz-plot anal- 
ysis will definitely keep us company over the next few 
years. There will be a lot of work for both experimen- 
talists and theorists alike: synergy will be invaluable. 
Some complications have already emerged, especially 
in the charm field, others, unexpected, will only be- 
come clearer when we delve deeper into the beauty 
sector. B s will be a completely new chapter. The 
analysis is challenging but there are no shortcuts to- 
ward ambitious and high-precision studies and, ulti- 
mately, to New Physics searches. 




1.25 L! 

(a) W(Ktt) (GeV/c ! 



(b) M(K7T) (GeV/c ! ) 



Figure 9: Comparison between the LASS 1=1/2 phase + 
ChPT (continous line) and the F-vector phases (with ±1 a 
statical error bars); a) total F-vector phase; b) I = 1/2 F- 
vector phase. The vertical dashed line shows the location 
of the Kn'. 



implied by the simple isobar fits. We have thus shown 
that whatever n is revealed by our D + — > K~tt + tt + 
results, it is the same as that found in scattering data. 
Consequently, our analysis supports the conclusions of 
[18] and [19]. 



4. Conclusions 

Dalitz-plot analysis represents a unique, powerful 
and promising tool for studying the Heavy Flavor de- 
cay dynamics. There is a recent, vigorous effort to 
perform amplitude analysis: a more robust formalism 
has been implemented, many channels have been in- 
vestigated. The beauty community can benefit from 
charm experience and expertise. The high statistic 
D+ -> K-TT+TT+ from FOCUS showed us that D- 
decay can also teach us about Kir interaction much 
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